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Abstract 

^^ ■ In this paper, we show our discovery that state-transition networks in several chaotic dynamical 

systems are "scale-free networks," with a technique to understand a dynamical system as a whole, 
which we call the analysis for "Discretized-State Transition" (DST) networks; This scale-free nature 
is found universally in the logistic map, the sine map, the cubic map, the general symmetric map, 

Q . the sine-circle map, the Gaussian map, and the delayed logistic map. Our findings prove that there 

is a hidden order in chaos, which has not detected yet. Furthermore, we anticipate that our study 
opens up a new way to a "network analysis approach to dynamical systems" for understanding 
complex phenomena. 
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I. INTRODUCTION 

Since time immeniorial, human beings have sought to understand the relationship be- 
tween "order" and "chaos." Some of the earhest conceptions associated order with pattern 
and consistency, and chaos with formlessness and fluidity. Scientific inquiries over the past 
four decades, however, have revealed that a number of nonlinear systems often exhibit ir- 
regular and aperiodic behaviors; their sequence of the states never repeat even though they 
are governed by a simple, deterministic rule [l|, |2|. This astonishing discovery of "determin- 
istic chaos" blurs the line between traditional conceptions of pattern and constancy, and 
formlessness and fiuidity. 

On one hand, a number of phenomena that behave chaotically have been found empirically 



in various systems in the natural, technological, and social world [3|-[l2j. On the other hand, a 
number of mathematical models exhibiting chaos have been studied, and consequently some 
identifiable orders have been recognized inside chaos from a certain viewpoint. A chaotic 
attractor visualization appears complex but coherent, for example the visually beautiful 
trajectories in a "phase space," which is a common representation for the evolution of a 



dynamical systems 
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14| . Besides, universal order has been found in bifurcation, which 



occurs when the control parameter of the system is changed [l5|, [l6| . These discoveries of 
order within chaos, however, focus only on a partial trajectory, namely attractor, rather 
than the whole behavior of systems. 

In this paper, we introduce a new technique that we call the analysis for "Discretized- 
State Transition" (DST) networks, in order to understand the whole behavior of dynamical 
systems rather than a partial behavior in an attractor. Then, with using the technique, we 
show that the state-transition networks in several chaotic dynamical systems are universally 



■'scale-free networks" 
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181], where there are quite many nodes having few links but also a 



few hubs that have an extraordinarily large number of links. 

II. DISCRETIZED-STATE-TRANSITION (DST) NETWORKS 

The technique for understanding the whole structure of behaviors in a dynamical system, 
the analysis for "Discretized-State Transition" (DST) networks, follows a two-step proce- 
dure: First, a continuous unit interval of variables is subdivided into a finite number of 
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FIG. 1. Building Discretized-State Transition (DST) networks: To obtain the DST networks for a 
dynamical system, the continuous unit interval of variables is subdivided into a finite number of 
subintervals, which is called "states," and then a directed network is build by connecting from one 
state to a successive state along the direction of the transition. 

subintervals, "states"; Second, a directed network is built with using these states, where 
a node represents a state and a link represents a transition from one state to a successive 
state. 

Concretely, first, we subdivide a unit interval of variables into a finite number of subin- 
tervals and call them "discretized states," or simply "states." The discretization is carried 
out by rounding the value of variables to d decimal places. Accordingly, one obtains a finite 
set of states whose size is given by A = 10"^. Mathematically, the map is described as 
a composite function as Xn+i = h{f{xn)) = h o f{xn), where the function h{x) is for the 
rounding function, which is either the round-up, round-off, or round-down function. 

Second, we build a directed network by connecting a link from one state to a successive 
state determined by the dynamical system. The order of the network A^, namely the number 



an initial condition 




FIG. 2. DST networks and Numerical Simulation. Numerical simulation shows an instance of flow 
on the DST networks 

of nodes, is equal to 1/A + 1. The out-degree of every node must be always equal to 1 since 
it is a deterministic system, while the in-degree of each node can exhibit various values. The 
whole behavior is often mapped into more than one connected components, which represent 
"basins of attraction." Each connected component must have only one loop or cycle, which 
represents "fixed point" or "periodic cycle," where either is an "attractor". Fig. [T] shows the 
simple schematic that summaries the technique, specifying the correspondent terms in the 
science of complexity with bracket under the terms in the network science. 

As a kind of directed network is occasionally called "river network" , it may be helpful 
to understand a state-transition network with a metaphor of "river." A node represents a 
geographical point in the river, and a link represents a connection from a point to another. 
The direction of a flow in the river is flxed. While there are no branches of the flow, however 
you can flnd confluences of two or more tributaries everywhere. To be exact, the network is 
not a usual river but "dried-up" river, where there are no flows on the riverbed. 

Determining a starting point and discharging water, you will see that the water flow 
downstream on the river network. That is a happening you witness when conducting a 
numerical simulation of system's evolution in time (Fig. [2]). Thus exploring a dynamical 
system with a numerical simulation means to observe an instance of flow starting from a 
point of the river, and then typical techniques of drawing a trajectory on phase space and 
bifurcation diagram are just to sketch out the instances. Our technique of a DST network, 
on the contrary, is to draw a "map" of the whole structure of the dried-up river networks, so 
one can overview the landscape from bird's-eye view. Thus the technique is much different 
from the typical technique of numerical simulations for understanding dynamical systems. 

Even though there have been some pioneering studies of dynamical systems with flnite- 
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FIG. 3. The DST networks for the logistic map with // = 1 in the case of round-up into d decimal 
places ranging from d = 1 to d = 4. 
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23| and studies of a state-transition network of discrete dynamical systems [2; 



29| , little is known about the nature of the state-transition networks in chaotic dynamical 



systems. 

III. SCALE-FREE DST NETWORKS IN THE LOGISTIC MAP 

We now move on to our discovery about chaotic dynamical systems with using the tech- 
nique; that is, state-transition networks in several chaotic dynamical systems are "scale-free 
networks." Such systems we found are mainly one-dimensional dynamical systems such as 
the logistic map, the sine map, the cubic maps, the general symmetric map, the Gaussian 
map, the sine-circle map, but also two-dimensional systems like the delayed logistic map. 

To demonstrate our discovery, we start with a logistic map, which may be the simplest 
and most well known map. The function for the logistic map is given by f{x) = 4/ix(l — x) 
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FIG. 4. The cumulative in-degree distributions of DST networks for the logistic map with /i = 1 
in the case from d = 1 to d = 7. The dashed line has slope 2, implying that the networks are 
scale- free networks with the degree exponent 7 = 1, regardless of the value of d. 



, where x is a variable between and 1, and yU is a control parameter ranging from to 1. 
Using this function, the value Xn+i is determined by the previous value Xn as Xn+i = fixn), 
and this is known as a simple population model with non-overlapping generations |5|. In 
spite of the simple looking of the function, the model exhibits various kinds of behavior 
according to the control parameter /i; for < /i < 0.75 , all initial conditions attract to the 
fixed point; for 0.75 < jj, < 0.89..., x oscillates as a periodic cycle; for 0.89... < /i < 1.00, the 
system exhibits chaos, namely the period is infinitely long. 

Fig. [3] shows DST networks of the logistic map with yU = 1 in the case of round-up into 
d decimal places, ranging from rf = 1 to d = 4, and fig. H] shows the cumulative in-degree 
distribution for the states whose in-degree /;; > , with d ranging from 1 to 7. The dashed 
line in the figure has slope 2, therefore distributions follow a power law with the exponent 
7 = 1 except in the regions where in-degree k is extremely high. Consequently, the result 
implies that these networks are "scale-free networks." 

Next, we investigated whether the scale- free nature of a DST network can be observed 
when varying the control parameter /x, which governs the behavior of the system. Fig. 
shows the cumulative in-degree distribution for DST networks of the logistic map with /^ 
ranging from to 1 by 0.125 in the case d = 7. For all values, the distributions follow a 
power law with the exponent 7 = 1 except in the regions where in-degree k is extremely 
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FIG. 5. The cumulative in-degree distributions of DST networks for the logistic map from /x = 0.125 
to 1.000 by 0.125 in the case d = 7. The dashed line has slope 2, implying that the networks are 
scale- free networks with the degree exponent 7 = 1, regardless of the value of the parameter fi. 

high. The result means that the scale-free nature is independent of the parameter value. It 
also implies that the scale-free nature of state-transition networks is ascribable to the map 
function, not the chaotic regime of the system's behavior. Fig. [6] shows DST networks in 
the case d = 3, which displays that the network structures are changed as the parameter n 
is varied, even though the scale-free property remains. 

Consequently, the results of our numerical explorations demonstrate that the DST net- 
works of the logistic map are scale-free networks, independent of the control parameter 
/i. 

IV. UNDERLYING MECHANISM FOR SCALE-FREE DST NETWORKS 



We turn now to examine the underlying mechanism with mathematical consideration. 
Investigating the correlation between the value of x and its in-degree k^ by numerical com- 
putation, it turns out that in-degree kr^ becomes higher as the value of x is getting to be 
larger. From the observation, it turns out that in-degree depends on the slope of the map 
function, and the hub states are located at the top of the parabola, whose slope is quite 
small. We shall describe a relation between y and x with using the size of subinterval A and 
the coefficient Cy, as shown in fig. [71 Based on the observation that the in-degree becomes 
higher as y is getting to be higher, the rank of state y is expressed as Vy = {fi — y)/A + 1. 
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FIG. 6. Structural change of the DST networks for the fogistic map in the case d = 3 with 
fi = 0.125, n = 0.250, /i = 0.375, fi = 0.500, /i = 0.625, /i = 0.750, ^ = 0.875, and /u = 1.000. 



In addition, as pointed out in some literatures [30|, ISlJ , cumulative degree distribution can 
be calculated from the rank / frequency plots as P^{k) = r/N, where r is the rank and N 
is a total number of nodes. 

Replacing r and A^ of the cumulative distribution function P^{k), after some mathe- 
matical deduction (see more details in Appendix |A]), we obtain the equations for describing 
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FIG. 7. The plot of the logistic map function with /x = 1. Note that the function is symmetric with 
respect to the critical point x* = 0.5, and then x„+i always yields two previous values of x„ but 
the case x„+i = 1. The box shows relation between a subinterval on y-axis and its corresponding 
range on the x-axis for the logistic map function. 

cumulative in-degree distribution of DST networks for the logistic map as follows: 



P>{k) oc {- + fiAk 



where in-degree k is in the range < k < k"^""^ and k"^"'^ = 1/ \/JiK. Recall that < ^A < 1, 
therefore the second term makes a large effect on the distribution only when k is quite 
close to k"^""^ . It means that the "hub" states will have higher in-degree than the typical 
scale-free network whose in-degree distribution follows a strict power law. Fig. [8] shows the 
fitting between the results of numerical computation and its mathematical predictions. As 
a result, we come to a conclusion that the underlying mechanism presented above is valid 
for the model to explain how scale-free networks are formed. 

Here, taking the limit of the equation, we have: 



limP>(A;) oc A;-^ 



The result shows that the cumulative in-degree distribution becomes a power-law distribu- 
tion in a strict sense, as the size of subintervals. A, is getting to be infinitely small. 
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FIG. 8. Comparison between the results of numerical computations and mathematical predictions 
with /i = 0.1 and /i = 1.0 in the case d = 6. 

V. SCALE-FREE DST NETWORKS IN OTHER CHAOTIC DYNAMICAL SYS- 
TEMS 



Since we focused only on the logistic map so far, we now broaden the scope to other 
maps in order to show that the scale-free DST networks are universally found in several 
chaotic dynamical systems. Fig. [9] shows the plot of the following map functions and their 
cumulative in-degree distributions: the sine map /(x) = fisin{Trx), and cubic maps /(x) = 
3\/3/xx(l —x'^)/2 and /(x) = 27/ix^(l — x)/4. The result implies that the scale-free property 
of DST networks is not specific to the logistic map. Furthermore, the similarity in the degree 
distributions and its exponents among several systems suggests a new kind of "universality" 
of state transitions, which is different from the well-known universality of bifurcation 15l.ll6|. 
Examples of the DST networks for these maps are shown in fig. [TOl 

Next, in order to understand how the flatness around the top of the parabola in a map 
function, we investigated the general symmetric map given by /(x) = fill — |2x — 1|"). 
Fig. [TT] shows the plot of the map function and their cumulative in-degree distributions, 
varying the value of a from 1.5 to 4.0 by 0.5. Note that the parameter /i is controlling the 
flatness around the top and the case /i = 2 is the same as the logistic map. The figure of the 
cumulative in-degree distribution shows that the parameter a influences the degree exponent 
7; however the scale-free property is maintained in the all case of a > 1.5. Examples of the 



10 



DST networks for the map are shown in fig. 

Finally, we explored the DST networks for the other types of map functions: exponential 
function, discontinuous function, and two-dimensional function (Fig. US]). First, Gaussian 
map given by /(x) = 6 + e~"^ , which is an exponential function, is scale-free network 



in the certain range of the parameter a. Second, Sine-circle map 32| given by f{x 



x + h— (a/2'7r)x {mod 1), which is discontinuous function, is also scale-free network. Third, 



the delayed logistic map 33| given by f{x,y) = (ax{l — y),x), which is two-dimensional 
function, is scale-free network as well. In all cases, the in-degree exponent is universally 
7 = 1. Examples of the DST networks for these maps are shown in fig. [TH 

VI. CONCLUSION 

Our discovery that state-transition networks of chaotic dynamical systems are organized 
in a certain way implies that our technique of mapping the transitions among discretized 



states is useful to capture t 
the network science 
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34 



le nature of the system as a whole. The recent development of 



36l | enables us to understand complex state-transition networks 
in light of a "scale-free" concept rather than random viewpoint. 

On the other hand, for network science, our study provides a demonstration of analyzing 
a network embedded in temporal behavior, where node is a "state" rather than either a sub- 
stance or an actor. Note that, in such a network, more than one node cannot exist together 
at the same time. Network science has not studied this kind of temporally-embedded net- 
works in much detail, and thus our study can be considered a step towards a new application 
of network analysis methods. 

Thus we anticipate that our study opens up a new way to a "network- analysis approach 
to dynamical systems," where one can understand complex systems with methods and tools 
developed in the network science. 
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FIG. 9. (a) The plots for the logistic map, the sine map, and the cubic maps, jjl and (12; and (b) 
their cumulative in-degree distributions in the case d = 6. The dashed line has slope 2, implying 
that the networks are scale- free networks with the degree exponent 7 = 1. 
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FIG. 10. Scale-free DST networks, in the case d = 3, for (a) the logistic map with fi = 1.0; (b) the 
sine map with /i = 1.0 ; (c) the cubic map jjl with /i = 1.0; and (d) the cubic map '^2 with /i = 1.0. 
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FIG. 11. (a) The plots for the general symmetric map with a ranging from 1.5 to 4.0 by 0.5; and 
(b) their cumulative in-degree distributions in the case d = 6. This figure shows that the parameter 
a influences the degree exponent 7, while the scale-free property is maintained. 









FIG. 12. Scale-free DST networks for the general symmetric map in the case d = 3 with (a) 
a = 1.5; (b) a = 2.0; (c) a = 2.5; (d) a = 3.0; (e) a = 3.5; and (f) a = 4.0. 
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FIG. 13. (a) The plots for the Gaussian map and the sine-circle map; and (b) their cumulative 
in-degree distributions for the Gaussian map and the sine-circle map in the case d = Q and the 
delayed logistic map in the case d = 3. The dashed line has slope 2, implying that the networks 
are scale- free networks with the degree exponent 7 = 1. 
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FIG. 14. Scale-free DST networks for (a) the sine-circle map with a = 4.0 and b = 0.5 in the case 
d = 3; (b) the Gaussian map with a = 1.0 and b = —0.3 in the case d = 3; and (c) the delayed 
logistic map with a = 2.27 in the case d = 2. 
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Appendix A: Mathematical Proof of Scale- Free DST Networks for the Logistic Map 

As shown in the paper, our simulation results demonstrate that the DST networks of the 
logistic map are scale-free networks with any values of control parameter, fi, and resolution 
for the discretization, A. Here we provide a mathematical proof of our findings. For the 
purpose, we shall begin with calculating the number of the states whose in-degree is higher 
than zero, and then calculating the cumulative in-degree distribution. 

1. Number of the states whose in-degree is higher than zero 

Suppose that / is the unimodal map of the interval [0, 1] that is concave downward, 
symmetrical with respect to the critical point, and /(O) = /(I) = (see Fig. [7]). While 
there are two solutions for the inverse function f~^, we shall indicate a solution on the left 
hand of parabola as f£^. Notice that, due to the symmetrical nature, in-degree k is always 
an even number but the critical point. 

In order to obtain the number of nodes who have incoming link, A^fc>0) we shall describe 
a relation between y and x with using the size of subinterval A and the coefficient Cy (see 
Fig. Wj. We then have y = f{x) and y — A = f{x — CyA), and, with the inverse function, 
X = fi^iy) and X — CyA = f]^^{y — A). Cancelling y from these equations gives the equation: 

_ x-fz\f{x)-A) 

y A 

On the other hand, cancelling x gives the equation: 

_ fl\y) - fl\y - A) 
Cy- ^ 

Accordingly, since ky = 2cy by definition, in-degree of the state y is calculated as: 

ky = l{fE'-fE\y-^))- 

Let the coordinates of the point where k = 4, namely Cy = 2, be (X, Y), and therefore it 
must satisfy that Y = f{X) and cy = (fE^ — fE^O^ ~ ^)) /^ = 2. We now consider the 
number of states for two ranges: Cy > Cy and < Cy < Cy. First, the number of states in 
the range Cy > cy is calculated as: 

15 



/x-r 



N, 



A 



Cy>CY 



A 



+ 1 



+ 1 if Y >0 



if Y<0. 



Second, the number of states in the range < Cy < cy is calculated as: 

X 



N< 



0<Cy<CY 



A 




+ 1 if X>0 
if X<0. 



Note that, since the map function is concave downward, /(O) = /(I) = 0, and focusing the 
left hand of the parabola, there are only two combinations of the conditions: X > A F > 
and X < A y < 0. Therefore, the number of states in the range < c, namely the number 
of the states whose in-degree is higher than zero, Nk>o, is given as follows: 

r fi + X-Y 

if X <0 A Y <0. 



2 i/X>OAr>0 



R. 

A 



2. Cumulative in-degree distribution 



As shown in some literatures |30|, l3l|, cumulative degree distribution can be calculated 
from the rank / frequency plots as: 

r 
P>(k) = — , 

where r is the rank and X is a total number of nodes. Based on the observation that the 
in-degree becomes higher as y is getting to be higher, the rank of state y is expressed as: 



A 



+ 1. 



Replacing r and N of the cumulative distribution function P^{k), we obtain: 

/i + A-y 



P>{k, 



/i + 2A + X - r 
/i + A-y 
I /^ + A 



i/X>OAF>0 

i/ x<o A r<o. 
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3. Cumulative in-degree distribution in the case of the logistic map 

In the case of the logistic map f{x) = 4yua;(l — x), a solution of the inverse function 
is given by f£^{y) = ( 1 — Jl — y/fi j /2, and the coordinates of the point where k = A, 
(X, Y), is given by: 



{X,Y) 



2 



1 + 16;uA| A 1 



16/x ' 2 64/i 

Solving the equation of ky for y, we obtain: 



fx - 4/iA^ 



1 



A 



y = ^ 



-2 + /iAA;2 . 



4/iA;2 4 

Since < A < 1/10 is always satisfied in the condition of this paper, the condition 
X > A F > means 1/(8 + 16A) < yU < 1; and the condition X < A F < means 
< yU < 1/(8 + 16A). Thus, the number of states whose in-degree is higher than zero is 
given by: 

' 1 - 4| - 1 + 16yuA| + 32/i(l + A + 8yuA2) 



X,, 



64yuA 



^/ 



1 
8+16A 



I A 



if < jj, < 



< /i< 1 



1 



16A 



Note that the maximum in-degree is observed at the highest point of the parabola y* = n, 
namely at the critical point x* = 1/2, and the maximum in-degree is estimated by: 



k' 



/JJA 



Consequently, in the case 1/(8 -|- 16 A) < fx < I, the cumulative in-degree distribution 
P^{ky) is given by: 



P>{k) 
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1 - 4| - 1 + 16/xA| + 32/i(l + 3A + 8/iA2) VA; 
In the case < /i < 1/(8 + 16A), the distribution is given by: 



- + fiAk 



P>{k) 



An{n + A) \k 



1 \ 2 1 - V A/yU + 1 
+ /iAA; if ^—:^ <k 



A 



if 0<k< 



1- JA//i + l 
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Thus, it is proved that the cumulative in-degree distribution of the DST network for the 
logistic map follows a law given by: 

P>(A;)oc f^ + fiAk) , 

where in-degree k is in the range < k < k"^""^ and /c™"^ = l/^/JxA. Recall that < fiA < 1, 
therefore the second term makes a large effect on the distribution, only when k is quite close 
to A;™-"^. It means that the "hub" states will have higher in-degree than the typical scale- 
free network whose in-degree distribution follows a strict power law. Fig. [8] shows the fitness 
between the results of numerical computation and mathematical predictions. 

In order to understand the distribution when the size of subintervals. A, is enough small, 
we take the limit of A into 0. In the case 1/(8 + 16A) < /i < 1, 

Ifi 

lim P>(k) = -k-^, 

A^o ^ ' -5 + 32/i + 128/i2 



and in the case < fj, < 1/(8 + 16A), 



lim P>(k) = < 



-^k-^ if —<k 

1 if < A;< — . 

■' - 2/i 



Thus, it is clear that the cumulative in-degree distribution of the DST network for the 
logistic map follows a strict power law when A is enough small as follows: 

P>{k) oc A;-l 

In conclusion, the DST networks for the logistic map become "perfect" scale-free networks 
whose degree distribution strictly follows a power law with the exponent 7 = 1, as the size 
of subintervals is getting smaller. 
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